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ABSTRACT 

We present observations of CLASS B2108+213, the widest separation gravitational 
lens system discovered by the Cosmic Lens All- Sky Survey. Radio imaging using the 
VLA at 8.46 GHz and MERLIN at 5 GHz shows two compact components separated by 
4.56 arcsec with a faint third component in between which we believe is emission from a 
lensing galaxy. 5-GHz VLBA observations reveal milliarcsecond-scale structure in the 
two lensed images that is consistent with gravitational lensing. Optical emission from 
the two lensed images and two lensing galaxies within the Einstein radius is detected in 
Hubble Space Telescope imaging. Furthermore, an optical gravitational arc, associated 
with the strongest lensed component, has been detected. Surrounding the system are 
a number of faint galaxies which may help explain the wide image separation. A 
plausible mass distribution model for CLASS B2108+213 is also presented. 

Key words: quasars: individual: CLASS B2108+213, cosmology: observations, grav- 
itational lensing. 



1 INTRODUCTION 

The Cosmic Lens All-Sk y Survey (CLASS; iMvers et afl 
120031: iBrowne et aD I2003T) is the largest statistically com- 
plete gravitational lens survey ever undertaken. A complete 
sample of 11 685 flat-spectrum radio sources (S v oc u a where 
a ^ —0.5 between 1.4 and 4.85 GHz) was observed with the 
Very Large Array (VLA) in A configuration at 8.46 GHz 
between 1994 and 1999. The 0.2-arcsec angular resolution 
provided by the VLA allowed the straight forward identi- 
fication of galaxy-scale gravitational lenses, which typically 
have image separations of ~ 1 arcsec. CLASS sources were 
identified as promising gravitational lens candidates if the 
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VLA 8.46 GHz image contained multiple compact compo- 
nents with Gaussian FWHM ^ 170 mas, the image sepa- 
rations were between 0.3 ^ Ad ^ 15 arcsec, the total inte- 
grated flux-density was 58.46 ^ 20 mjy and the component 
flux-density ratio was 10 : 1. The detection of compact 
structure and extended emission in the multiple components 
at higher resolutions with MERLIN (Multi-Element Radio 
Linked Interferometer Network) and the VLBA (Very Long 
Baseline Array) confirmed the lensing hypothesis for 21 
gra vitational lenses dur ing the first three phases of CLASS 
(see IBrowne et al"ll2003l and references therein for a discus- 
sion of all the gravitational lenses discovered by CLASS). 

In this paper, we present observations of the gravita- 
tional lens candidate CLASS B2108+213, which was discov- 
ered during the fourth and final phase of CLASS. In Section 
1^1 we describe the VLA, MERLIN and VLBA radio observa- 
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tions which led to the discovery of the gravitational lens can- 
didate CLASS B2108+213. Follow-up optical and infrared 
imaging of the system and the surrounding field, taken with 
the Hubble Space Telescope, is presented in Section Us- 
ing the available observational data, the lensing hypothesis, 
along with an analysis of the lensing potential, is discussed 
in Section |3] Finally, in Section |5] we summarise our find- 
ings and outline the future work we intend to carry out on 
this intriguing new gravitational lens from CLASS. For all 
calculations we adopt an Q.m = 0.3, = 0.7 flat-universe, 
with a Hubble parameter of Ho = 70 km s -1 Mpc -1 . 



2 RADIO OBSERVATIONS 
2.1 VLA 8.46 GHz observation 

CLASS B2108+213 was observed with the VLA in A con- 
figuration at 8.46 GHz as part of the final phase of CLASS 
on 1999 August 17. The total on-source integration time 
was ~ 30 seconds. A 3.3 second correlator integration time 
with two 50 MHz IFs were used. 3C286 was chosen as the 
primary flux-density calibrator and phase-referencing was 
carried out with J VAS B2059+034 (Jodrell Bank VLA As- 
trometric Survey ; [Patn aik et al1ll992t iBrowne et al J TT998t 
IWilkinson et al J |j99^L The data were edited and cali- 
brated in the standard way using AIPS (Astronomical image 
processing Software) and imaged wit hin the Caltech VLBI 
difference mapping package (DlFMAP; |SherAe rdlll997h using 
the CLASS mapping script jMvers etall20oljfr TIie itu-data 
were naturally weighted and elliptical Gaussian model com- 
ponents were fitted to each image. The rms noise level was 
180 fiJy beam -1 and the beam size was 0.26 x 0.23 arcsec 2 . 

The VLA detected three compact components, shown in 
Fig. Q whose flux-densities and relative positions are given 
in Table and Table |5J respectively. With an image sepa- 
ration between components A and B of 4.56 arcsec, CLASS 
B2108+213 has a wider image separation than any grav- 
itational lens system thus far discovered by CLASS (see 
Fig. [5] for the CLASS gravitational lens image separation 
distribution). The total integrated flux-density of the sys- 
tem at 8.46 GHz is 18.5 mjy and the flux-density ratio be- 
tween the brightest and faintest components is ~ 8 : 1. 
Therefore, CLASS B2108+213 failed to meet the strict 
lens candidate selection criteria due to its total integrated 
8.46 GHz flux-density being less than 20 mjy. However, 
CLASS B2108+213 was followed up as a promising lens can- 
didate even though it was no longer part of the statistically 
well defined CLASS sample. 



VLA 8.46 GHz 17 August 1999 
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Figure 1. VLA 8.46 GHz observation of CLASS B2108+213 
taken during the final phase of CLASS on 1999 August 17. 
Three compact components were detected. The separation be- 
tween components A and B is 4.56 arcsec. Contours are at 
(-3,3,6,12,24,48) X 180 /uJy beam" 1 . The beam size is 259 X 
233 mas 2 at a position angle of 50.6°. North is up and east is left. 
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Figure 2. The image separation distribution of the 22 CLASS 
gravitational lenses. CLASS B2108+213, with a maximum image 
separation of 4.56 arcsec, is the widest separation gravitational 
lens system in the CLASS sample. 



2.2 MERLIN 5 GHz observations 

In order to determine if any of the compact components 
detected by the VLA showed extended emission at higher 
resolutions we observed CLASS B2108+213 with MER- 
LIN at 5 GHz for ~ 1.5 hours on 2000 March 10 and 
~ 9 hours on 2001 August 10. The same set-up was used 
for both observations. 3C286 was used as the primary flux- 
density calibrator and JVAS B2103+213 was used for phase- 
referencing. The point source OQ208 was used to deter- 
mine the antenna-based phase corrections. The observations 



were carried out using the left and right circular polariza- 
tions with 16 x 1 MHz channels. The data were initially 
edited and flux-density calibrated using the standard MER- 
LIN data reduction programmes dplot and tdproc before 
being processe d within AIPS usin g the MERLIN automated 
pipeline jThomasson et al"lll994T) . As before, the mapping 
was carried out within difmap, using natural weighting. The 
rms noise level and beam size were 298 /xJy beam -1 and 
63 x 49 mas 2 on 2000 March 10, and 109 ^Jy beam -1 and 
70 x 44 mas 2 on 2001 August 10. The observing conditions 
on 2001 August 10 were poor. As a consequence, an uncer- 
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Date 


Instrument 


Frequency 


Integration 


Beam size 




Flux density 










[GHz] 


time [min] 


[arcsec 2 , P. A.] 


A [mjy] 


B [mjy] 


C [mjy] 


[fijy beam ] 


1999 Aug 1 / 


VLA 


8.460 


0.5 


0.259 X 0.233, +51 


11.7 + 0.6 


5.3 + 0.3 


1.5 + 0.2 


180 


2000 Mar 10 


IVl-Llv LI IN 


4.994 


100.0 


0.0o3 X 0.049, +32 


8.8 ±0.5 


4.4 ± 0.4 


< 1.5 


298 


2001 Mar 21 


VLBA 


4.987 


100.0 


0.004 X 0.002, -18° 


6.8 ± 0.4 


2.6 + 0.2 


< 0.7 


130 


2001 Aug 10 


MERLIN 


4.994 


540.0 


0.070 x 0.044, +25° 


10.6 ± 1.6 


5.0 + 0.8 


1.1 + 0.2 


109 


2001 Doc 31 


VLA 


4.860 


2.0 


28.10 x 15.00, -66° 


30.0 + 1.5 






173 


2001 Dec 31 


VLA 


8.460 


2.0 


14.90 x 8.390, -60° 


25.8 ± 1.3 






110 


2001 Dec 31 


VLA 


14.940 


4.0 


9.500 x 4.810, -66° 


21.9 ± 1.1 






303 


2001 Dec 31 


VLA 


22.460 


7.0 


5.950 x 3.130, -61° 


10.6 ± 0.6 


7.3 + 0.5 




305 



Table 1. A summary of the radio observations of CLASS B2108+213. Non-detections are quoted at the 5<r level. For those observations 
where the resolution was insufficient to separate components the flux-density is listed under the brightest component and the others are 
marked <— . 



Date 


Instrument 


Comp. 


Act 


A<5 








[mas] 


[mas] 


1999 Aug 17 


VLA 


A 


+ 13 


+ 13 






B 


2139 ± 15 


-4027+ 11 






C 


1425 ± 34 


-2892 ± 34 


2001 Mar 21 


VLBA 


A 


0+1 


0+1 






B 


2131 ± 1 


-4032 ± 1 


2001 Aug 10 


MERLIN 


A 


+ 11 


0+11 






B 


2122 ± 11 


-4039 ± 11 






C 


1422 ± 13 


-2888 ± 13 



MERLIN 5 GHz 10 August 2001 



Table 2. The relative positions of each radio component of 
CLASS B2108+213 as measured with the VLA, VLBA and MER- 
LIN. The uncertainties in the relative positions have been calcu- 
lated from the uncertainty in the position of the Gaussian fitted 
to each component (i.e. beam size/SNR). The J2000 position of 
component A, as measured with the VLBA, is 21 h 10 ra 54. 0242 s 
+21°31'0.690". The flux-densities of each component can be 
found in Table Q] 



tainty of ~ 15 per cent has been assigned to the flux-density 
calibration for that observation. 

Components A and B were detected by the 2000 March 
10 observation and were found to be compact (Gaussian 
FWHM of 34 and 30 mas, respectively). However, com- 
ponent C was not detected. This was due to the radio 
emission from component C being extended and resolved 
out by the MERLIN beam, or having a flux-density be- 
low S4.994 < 1.5 mjy, the 5a image sensitivity. To inves- 
tigate these possibilities the deeper 5 GHz observation on 
2001 August 10 was undertaken. All three components of 
CLASS B2108+213 that were identified by the VLA were de- 
tected (see Fig. [3J| and fitted with single elliptical Gaussian 
model components. As previously observed, components A 
and B were found to have compact structures (Gaussian 
FWHM of 27 and 21 mas, respectively) consistent with lens- 
ing. However, component C was resolved (Gaussian FWHM 
of 64 mas). The flux-densities of each component from both 
observations are given in Tableland their relative positions, 
taken from the 2001 August 10 observation, are presented 
in Table H 



Relative R.A. (arcsec) 

Figure 3. MERLIN 5 GHz long-track observation of 
CLASS B2108+213 taken on 2001 August 10. The FWHM of the 
elliptical Gaussian model fitted to component C is comparable to 
the MERLIN beam size. Contours are at (-3, 3, 6, 12, 24, 48) X 
109 fijy beam -1 . The beam size is 70.1 X 44.1 mas 2 at a position 
angle of 24.6°. North is up and east is left. 



2.3 VLBA 5 GHz observation 

A snapshot observation with the VLBA at 5 GHz was made 
of CLASS B2108+213 on 2001 March 21. To provide uni- 
form leu-coverage the observation was split into 5 separate 20 
minute on-source integrations. Phase-referencing and fringe- 
fitting were carried out using JVAS B2103+213 and 3C345, 
respectively. The data were recorded in the left circular po- 
larization through 4 IFs, each with 8 MHz bandwidth. 1 Since 
the separation of components A and B was known to be 
4.56 arcsec we obtained two correlations for the data with 
each centred on the respective components. The calibration 
and data editing was performed in the standard way us- 
ing AlPS before being averaged into 2 MHz frequency chan- 



1 Note that each IF was split into 16 X 0.5 MHz channels during 
the data processing at the correlator. 
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nels and 10 second integrations. The data were cleaned and 
mapped in AIPS using imagr which allowed the data to be 
imaged without further frequency channel averaging. The 
rms noise level was 142 and 124 /xjy beam - for the respec- 
tive maps of components A and B, which are given in Fig. 
|U The beam size was 4x2 mas . 

As before, both A and B components were detected. 
The brighter component (A) contains weak extended emis- 
sion stretching to the south-west whereas component B is 
compact. Furthermore, the extension observed in compo- 
nent A is perpendicular to the axis connecting components 
A and B. This is consistent with the lensing hypothesis. 
However, component C was not detected. The sensitivity of 
these deep VLBA data should have been sufficient to detect 
component C if it were compact. The individual component 
flux-densities are given in Table Q and the relative positions 
are given in Table [5] 

2.4 VLA multi-frequency observations 

Flux-density measurements of CLASS B2108+213 at 4.86, 
8.46, 14.94 and 22.46 GHz were carried out with the VLA in 
D configuration on 2001 December 31. CLASS B2108+213 
was observed for 2 minutes at 4.86 and 8.46 GHz, 4 minutes 
at 14.94 GHz and 7 minutes at 22.46 GHz. A 3.3 second 
correlator integration time was used with two 50 MHz IFs. 
3C286 was used as the primary flux-density calibrator and 
phase-referencing was carried out using J2115+295, which 
was selected from the VLA calibrator list. The beam size 
and rms map noise at each observing frequency (4.86, 8.46, 
14.94 and 22.46 GHz) were 28.1 x 15.0, 14.9 x 8.4, 9.5 x 
4.8, 6.0 x 3.1 arcsec 2 and 173, 110, 303, 305 /LtJy beam" 1 , 
respectively. 

The data were calibrated and edited in the standard 
way using AIPS. Mapping and self-calibration were carried 
out within difmap. Natural weighting was used to maxi- 
mize the signal-to-noise of the maps. For the 4.86, 8.46 and 
14.94 GHz observations CLASS B2108+213 was unresolved, 
therefore a single elliptical Gaussian model component was 
fitted to the data. However, the beam size of the 22.46 GHz 
observation was sufficient to resolve component A from com- 
ponents B and C. Therefore, two elliptical Gaussian model 
components were fitted to the positions of A and B/C. The 
angular resolution was not sufficient to resolve components 
B and C. At each observing frequency, strong emission was 
detected from CLASS B2108+213 (see Table □ for the re- 
spective flux-densities) . 



3 OPTICAL & INFRARED OBSERVATIONS 

An important step in the lens candidate confirmation pro- 
cess is the optical identification of a suitable lensing galaxy 
and the lensed images. This is especially important for the 
case of CLASS B2108+213 since the wide image separation 
of 4.56 arcsec probably requires the presence of a group or 
cluster of lensing galaxies. 

High resolution imaging of CLASS B2108+213 and the 
surrounding field was carried out with the Hubble Space 
Telescope (GO-9744) on 2003 September 09. The optical 
data were taken with the Advanced Camera for Surveys 
(ACS), in the Wide Field Channel (WFC) mode, which has a 



field of view of 202 x 202 arcsec 2 and a pixel scale of 0.05 arc- 
sec pixel -1 . The total on-source integrations, taken through 
the F555W and F814W filters, were 2.2 ks and 2.3 ks, re- 
spectively. Infrared imaging through F160W was carried out 
with the Near Infrared Camera/Multi-Object Spectrograph 
(NICMOS). This 2.6 ks observation was taken with the NIC2 
camera in the multi-accum data taking mode. NIC2 has 
a field of view of 19.2 x 19.2 arcsec 2 and a pixel scale of 
0.075 arcsec pixel - . 

The data were reduced using standard stsdas tasks 
within iraf. 2 The ACS data for each filter were first 
processed using calacs, before being combined and cor- 
rected for the ACS geo metric distortion using multidrizzle 
jKoekemoer et alJl20oA . The infrared NICMOS data were 
reduced using the calnica and calnicb tasks. Further im- 
age processing (i.e. bad pixel and residual cosmic ray re- 
moval) was carried out using the fillbad and clean tasks 
within KAPPA and FIGARO, 3 respectively. A median weighted 
mosiac of the background-subtracted dithered images was 
created using ccdpack. The ACS F555W and F814W, and 
NICMOS F160W images of CLASS B2108+213 are pre- 
sented in Fig. |S] (top) . 

Optical and infrared emission coincident with the two 
lensed images (A and B) has been detected. Both lensed im- 
ages appear to be strong point sources, however, extended 
arc-like emission, which is presumably from the host galaxy, 
was also detected in the F814W and F160W images of com- 
ponent A. A lensing galaxy (Gl) has been found within 
6 mas of the position of radio component C (based on the 
MERLIN 5 GHz position). Faint emission from a fourth 
component (G2) within the Einstein radius of the system is 
likely to be from a companion lensing galaxy. Both deflectors 
have a morphology which is consistent with an early-type 
galaxy, although there is evidence that Gl may be disturbed 
(see below for surface brightness profile fitting). In Fig.||j]wc 
show the ACS F814W image of the CLASS B2108+213 field. 
Surrounding the lens system are a scatter of faint galaxies 
which may be contributing to the gravitational potential and 
might help account for the wide image separation. However, 
we defer any discussion of a possible group associated with 
CLASS B2108+213 to a follow-up paper (McKean et al., in 
preparation) . 

In order to carry out a photometric and astrometric 
analysis of the data, the surface brightness profiles of the 
two lensing galaxies were fitted with a de Vaucouleurs' pro- 
file, convolved with the ACS/NICMO S point spread func- 
tion, 4 using GALFIT iPeng et al.l 12002*). The parameters of 
the profiles fitted to the ACS F814W data-set, which has the 
highest signal-to-noise ratio of the three images, are given 
in Tableland the residual images, with the lensing galaxies 
subtracted, are shown in Fig. (bottom). The fitted profiles 
are a good representation of the galaxy cores, with relative 
residuals (i.e. observed/model) of ~ 4 percent. However, one 

2 IRAF (image Reduction and Analysis Facility) is distributed by 
the National Optical Astronomy Observatories, which are oper- 
ated by AURA, Inc., under cooperative agreement with the Na- 
tional Science Foundation. 

3 KAPPA, FIGARO and CCDPACK and are part of the STARLINK 
project - |http: //www.starlink.rl.ac.uk| 

4 The ACS and NICMOS point spread functions were generated 
using TINYTIM iKrisdll99^ . 
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Figure 4. Left: VLBA 5 GHz imaging of component A taken on 2001 March 21. Contours are at (-3,3,6,12,24) X 142 /ijy beam" 1 . 
The beam size is 4 X 2 mas 2 at a position angle of —18.1°. Right: The corresponding VLBA image of component B taken on 2001 March 
21. Contours are at (—3,3,6,12,24) X 124 /ijy beam" 1 . The beam size is 4 X 2 mas 2 at a position angle of —18°. North is up and east 
is left. 



of the most striking features of the lensing galaxy subtracted 
images is the asymmetry in the residuals in the outer tail of 
Gl. The fitted surface brightness profile over (under) repre- 
sents the observed profile toward the west (east). The likely 
cause of these asymmetric residuals is the interaction be- 
tween the lensing galaxies which could lead to the disturbed 
morphology of Gl. 

The magnitudes of the lensed components were mea- 
sured within elliptical apertures using the images with the 
lensing galaxies subtracted. The same size of aperture was 
used to calculate the magnitudes of A and B in the F555W, 
F814W and F160W images. The positions of the lensed com- 
ponents were found by fitting point sources to the strong 
core emission detected in each image. The pixel scales for 
each instrument, which are required for the astrometric anal- 
ysis of the data, were established by assuming that the ra- 
dio and optical core emission from components A and B 
was coincident. Using the precise AB image separation pro- 
vided by the VLBA 5 GHz observations, we found the (x, 
y) pixel scales (mas pixel -1 ) for the ACS and NICMOS to 
be (50.320, 49.983) and (75.929, 75.000), respectively. The 
relative positions and magnitudes of each component are 
presented in Table [3] Note that the magnitudes quoted in 
Table El have not been corrected for galactic extinction. 



CLASS B21084-213 FIELD HST/ACS F814W 
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Relative R.A. (arcsec) 
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Figure 6. Hubble Space Telescope ACS F814W imaging of the 
CLASS B2108+213 field. Surrounding the lens system are a scat- 
ter of faint galaxies which could be contributing to the lensing 
potential. North is up and east is left. 



4 DISCUSSION 

The observational evidence presented above will now be used 
to argue for the gravitational lensing hypothesis in CLASS 
B2108+213 and investigate the nature of the third radio 
component. A plausible model for the lensing potential is 
also presented. 



4.1 Evidence for the lensing hypothesis 

The radio data presented in Section [5] provide compelling 
evidence for the lensing hypothesis in CLASS B2108+213. 
Firstly, the VLA and MERLIN images have identified two 
compact components (A and B) separated by 4.56 arcsec 
with flat radio spectra, af' 46 ~ 0.15. High resolution imag- 
ing with the VLBA has shown that component A exhibits 
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Figure 5. Optical and infrared imaging of CLASS B2108+213 taken with the Hubble Space Telescope on 2003 September 09. Both 
lensed images (A and B) and two lensing galaxies (Gl and G2) have been detected. Top, from left to right, are the ACS F555W, ACS 
F814W and NICMOS F160W images, respectively. The same images, with the lensing galaxies subtracted, are shown on the bottom. The 
fitted galaxy profiles are a good approximation of the core emission from both lensing galaxies (i.e. observed/model relative residuals of 
~ 4 percent). In the ACS F555W image the lensed components appear point-like. Although a strong point source was detected, a faint 
underlying arc feature, which is presumably from the host galaxy, can also be seen in the ACS F814W and NICMOS F160W images of 
component A. North is up and east is left. 



Filter 


Comp. 


Aa 


AS 


Magnitude 






[mas] 


[mas] 


[Vega] 


F555W 


A 


Oil 


Oil 


22.56 ±0.15 




B 


2131 ±3 


-4032 ± 4 


23.73 ±0.17 




Gl 


1427 ± 5 


-2890 ± 5 


19.91 ±0.11 




G2 


1273 ± 5 


-1786 ± 5 


22.59 ±0.15 


F814W 


A 


Oil 


Oil 


20.96 ±0.15 




B 


2131 ± 2 


-4032 ± 3 


21.97 ±0.17 




Gl 


1428 ± 2 


-2888 ± 2 


17.63 ±0.11 




G2 


1275 ± 2 


-1786 ± 2 


20.61 ±0.15 


F160W 


A 


Oil 


Oil 


18.14 ±0.15 




B 


2131 ± 2 


-4032 ± 3 


18.86 ±0.17 




Gl 


1421 ± 2 


-2888 ± 2 


15.65 ±0.11 




G2 


1263 ± 1 


-1795 ±2 


18.36 ±0.15 



Table 3. The magnitudes (Vega) and relative positions of each 
optical component of CLASS B2108±213. 



extension stretching toward the south-west, which is perpen- 
dicular to the axis connecting components A and B; whereas 
the weaker component, B, is compact. This is typical of the 



Parameter Gl G2 



Total magnitude 


17.63 


20.61 


Effective radius 


1.36 arcsec 


0.19 arcsec 


Axial ratio 


0.86 


0.79 


Position angle 


57° 


134° 



Table 4. The model parameters of the de Vaucouleurs' pro- 
files fitted to the ACS F814W data of Gl and G2 using GALFIT 
JPeng et al . 2002). These parameters are consistent with those fit- 
ted to the ACS F555W and NICMOS F160W data. The position 
angles arc measured cast of north. 



gravitational lensing phenomena as surface brightness in the 
lense d images is conserved (e.g. ISchneider. Ehlers fc Falcol 
1992). Also, as gravitational lensing is achromatic and radio 
emission is not affected by dust absorption, we would expect 
the flux-density ratios of multiple images of the same back- 
ground source to be nearly identical over all frequencies. The 
flux-density ratios of components A and B at 5 and 8.46 GHz 
are S 5 , B /S 5 , A = 0.47±0.10 and S S ag,b / S 8 a6,a = 0.45±0.03 
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respectively. Furthermore, from the Hubble Space Telescope 
NICMOS F160W data presented in Section the infrared 
flux-density ratio of components A and B was found to be 
Sfwow,b I Sfwow,a — 0.51 ±0.04. Such similar flux-density 
ratios between 6 cm and 1.6 ^im provides strong confirma- 
tory evidence that components A and B are lensed images. 
Further evidence for the lensing hypothesis is provided by 
the under-lying gravitational arc detected in the F814W and 
F160W images of component A. The identification of two 
suitably positioned lensing galaxies also adds to the case for 
components A and B being images of the same background 



4.2 The nature of the third radio component 

One of the most interesting features of CLASS B2108+213 
is the possible detection of a third lensed image (although, 
the observations presented here suggest the emission from 
component C is probably from the lensing galaxy, Gl). 
Lensing theory predicts that for extended mass distribu- 
tions there should exist an odd (i.e. 3rd, 5th, etc) lensed 
image which i s positioned very close to the centr e of the lens- 
ing potential JPver fc Roederll98d : lBurkeil98 ll) . Therefore, 
the detection (or non-detection) of such an image can be 
use d to constrain the c ore-properties of early- type galax- 
ies jRusin fc Mai 1200 ll : iKeetorJ 120031) . Unfortunately, be- 
ing positioned near to the lens centre means that odd im- 
ages are also highly de-magnified. As a consequence their 
detection is extremely rare. The triple component gravi- 
tational lens system PMN J 1632-00 33 dWinn et alJl2002t 
IWinn. Rusin fc Kochanekl2003l.r2004t) currently has the best 
core image candidate. Other gravitational lenses which were 
found to have a candidate odd-image have subsequently 
shown that the emission is probably not from a lensed im- 
age but from a, lensi ng galaxy (e.g. IChen fc Hewitti Tl993t 
iFassnacht et al]ll999l) . 

The latter is almost certainly the case in CLASS 
B2108+213 for two reasons. Firstly, from the long-track 
MERLIN 5 GHz image presented in Section ^. 21 it is appar- 
ent that component C is extended (recall that the FWHM of 
the elliptical Gaussian fitted to C is comparable to the MER- 
LIN beam size) whereas components A and B are both com- 
pact. This conclusion is consistent with the non-detection of 
component C with the VLBA at 5 GHz. Secondly, the im- 
plied radio spectrum of component C differs from that of A 
and B. As can be seen in Fig. component C has an in- 
verted radio spectrum between 5 and 22.46 GHz, whereas A 
and B show a turnover. 6 Therefore, the component C radio 
emission is probably from the lensing galaxy, Gl, because 
the emission is not compact and has a different radio spec- 
trum from the lensed images, A and B. There is a subtle flaw 
to this argument. The implied radio spectrum of component 
C is rising at 22.46 GHz. This suggests that the source of 
the component C radio emission is very compact. However, 
the MERLIN 5 GHz imaging has shown component C to 



5 This flux-density ratio has been calculated using only the emis- 
sion from the strong point sources. 

6 The flux-density of component C at 22.46 GHz has been es- 
timated by assuming that the flux-density ratio of A and B at 
22.46 GHz is the same as at 5 and 8.46 GHz. 




10 100 
Frequency (GHz) 

Figure 7. The radio spectral energy distribution of CLASS 
B2108+213 and co mponents A, B and C . The em pty square and 
circle are the NVSS iCondon et alll99Sl) and GB6 iGregorv et alJ 
Il99tt flux-densities of CLASS B2108+213 respectively. The 
empty and filled triangles represent the CLASS VLA and multi- 
frequency VLA flux-densities. The flux-densities measured by 
MERLIN are the filled squares. Notice that the GB6 and VLA 
~ 4.85 GHz flux-densities are identical. The 22.46 GHz flux- 
densities for components B and C have been estimated by as- 
suming that the flux-density ratio, Sb/Sa is constant. Note that 
the sum of the radio component flux-densities is less than the to- 
tal flux-density measured by GB6 and the VLA multi-frequency 
observations. This is most apparent at ~ 4.85 GHz and is proba- 
bly due to resolution effects. 



be extended. A possible explanation for this contradiction 
is scattering in the lensing galaxy, which cou ld decrease the 
surface brightness of the radio emission (e.g. iGuira do et alJ 
ll999UBiggs et ail2003UWinn et aJ2003UBiggs et all2' 
However, accepting a scatter broadening explanation for the 
observed extension opens up the possibility that component 
C may be a lensed image since surface brightness will no 
longer be conserved. 

Further radio observations to determine the actual ra- 
dio spectra of the CLASS B2108+213 radio components will 
need to be carried out to resolve this issue. However, the 
most likely scenario is that the radio emission from compo- 
nent C is from the lensing galaxy. Although this conclusion 
rules out the exciting possibility of using CLASS B2108+213 
for studying the core properties of the lensing potential, an 
accurate position for the lensing galaxy Gf has probably 
been obtained. This will improve the mass model of CLASS 
B2108+213 which is discussed in the following section. 



4.3 Modelling the lensing potential 

Using the observational data presented in this paper, a pre- 
liminary analysis of the CLASS B2108+213 lensing potential 
can now be carried out. The lensed image positions and flux- 
densities provided by the VLBA 5 GHz data (see Tableland 
Table |5J give five model constraints; two from each image 
position and one from the flux-density ratio. The positions 
of the lensing galaxies have been taken from the NICMOS 
F160W data and allowed to vary within the quoted uncer- 
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Comp. Parameter 



Gl 


a, 8 


1421, -2888 mas 




bci 


1786 mas 


G2 


a, 8 


1263, -1795 mas 




bQ2 


518 mas 


Shear 


Magnitude 


0.018 




Position angle 


118.1° 


A 


a, 8 


0, mas 




Flux density (5 GHz) 


6.7 mjy 


B 


a, 8 


2131, -4032 mas 




Flux density (5 GHz) 


2.6 mjy 


Source 


a, 8 


1028, -2033 mas 




Flux density (5 GHz) 


1.4 mjy 



Table 5. The parameters of the CLASS B2108+213 mass model. 
The model successfully reproduces the observed positions and 
flux-densities of the lensed components (see Tables HI and l2t . All 
positions are measured relative to component A. 



tainties (see Table [3J. Each lensing galaxy is represented by 
a singular isothermal sphere, which has a convergence (k) 
of, 



where b is the Einstein radius and r is the radial distance 
from the lens centre. For a singular isothermal sphere the 
Einstein radius is given by, 




where a v is the one-dimensional velocity dispersion, and 
Dig and D s are the angular diameter distances to the back- 
ground source from the lens and observer, respectively. The 
ratio of the Gl and G2 luminosities (L) has been used to 
constrain the mass ratio of the lensing galaxies (&G2/&G1 = 
0.29 ± 0.05) by invoking the Fab er- Jackson relation (L o c b 2 
for a singular isothermal sphere; iFaber fc Jacksonlll976l) . A 
mass model containing two singular isothermal spheres, with 
the mass ratio fixed, requires three free parameters; two for 
the source position and one for the Einstein radius of Gl. 
The final two free parameters have been used for the position 
angle and strength of an external shear. 

The time-delay surface of our mass model for CLASS 
B2108+213 is shown in Fig. |S| and the model parameters 
are given in Table |S] We find that two singular isothermal 
spheres, with an external shear component, reproduces the 
observed positions and flux-densities of the lensed images. 

Due to the limited number of observational constraints, 
this mass model has no degrees of freedom. Further improve- 
ments of this simple model, which only serves to show that 
the system can be explained by gravitational lensing, will be 
presented in a companion paper (McKean et al., in prepara- 
tion). Additional constraints from stellar kinematic data for 
Gl, the extended structure of the lensed images, and redshift 
and luminosity data on putative group members around the 
lens system will be used to further investigate the CLASS 
B2108+213 lensing potential. 




Arcsec 

Figure 8. The time-delay surface for the CLASS B2108+213 
mass model. The model consists of two singular isothermal 
spheres and an external shear. The time-delay is measured in days 
and assumes that G2 is at the same redshift as Gl (zqi = 0.365; 
McKean et al., in preparation) and that the source redshift is 
1.5. However, since the source redshift is unknown, the cosmol- 
ogy dependent time-delay is only an estimate. For different source 
redshifts (or Hubble constant) , each contour level scales by a con- 
stant factor. The model parameters are given in Table 151 



5 CONCLUSIONS AND FURTHER WORK 

We have presented radio and optical observations of the new 
gravitational lens CLASS B2108+213. This system com- 
prises two lensed images, which are separated by 4.56 arcsec, 
thus making CLASS B2108+213 the widest separation grav- 
itational lens discovered by CLASS. Optical imaging found 
two lensing galaxies within the Einstein radius of the system. 
Intriguingly, we also detected radio emission from a third 
component which was coincident with the optical emission 
from the primary lensing galaxy. The extension of the third 
radio component and its different radio spectrum from the 
two lensed images suggests that it is unlikely to be a third 
lensed image. However, it is too early to completely dis- 
count the possibility that scatter broadening effects might 
have been responsible for modifying the properties of a third 
image. The positions and flux-densities of the lensed images 
were successfully reproduced by a model containing two sin- 
gular isothermal spheres and an external shear. However, 
further observational constraints will need to be obtained 
for a detailed analysis of the CLASS B2108+213 lensing po- 
tential to be carried out. 

The wide image separation of CLASS B2108+213 im- 
plies a group/cluster assisted lensing potential. Therefore, 
CLASS B2108+213 provides an excellent opportunity to 
study the structure and mass distribution of galaxy groups 
beyond the local universe (z > 0.2). We have begun a spec- 
troscopic programme to determine the redshifts of the lens- 
ing galaxies and lensed source in CLASS B2108+213, and 
to identify which nearby galaxies are part of the lensing 
group. Thus far we have found four nearby galaxies at simi- 
lar redshifts to Gl (zgi = 0.365; McKean et al., in prepara- 
tion). Furthermore, our spectrum of Gl has sufficient signal- 
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to-noise for the stellar velocity dispersion to be measured. 
These observations will provide additional constraints for a 
future mass model. 

The other interesting feature of CLASS B2108+213 is 
the possible detection of a third lensed image. Although the 
observations presented here show that this is almost cer- 
tainly not the case, this conclusion could be confirmed by 
measuring the flux-density, polarization and size of the radio 
emission from each component as a function of frequency. If 
scattering of a third lensed image is occurring then we would 
expect the high frequency spectrum of each radio component 
to be identical. 
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